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We study supersymmetric (SUSY) responses to a photoassociation (PA) process in a mixture 
of Bose molecules h and Fermi atoms / which turn to mutual superpartners for a set of proper 
parameters. We consider the molecule 6 to be a bound state of the atom / and another Fermi 
atom F with different species. The h-f mixture and a free F atom gas are loaded in an optical 
lattice. The SUSY nature of the mixture can be signaled in the response to a photon-induced 
atom-molecule transition: While two new types of fermionic excitations, an individual h particle-/ 
hole pair continuum and the Nambu-Goldstone-fermion-like ( or goldstino-like) collective mode, are 
concomitant for a generic b-f mixture, the former is completely suppressed in the SUSY b-f mixture 
and the zero-momentum mode of the latter approaches to an exact eigenstate. This SUSY response 
can be detected by means of the spectroscopy method, e.g., the PA spectrum which displays the 
molecular formation rate of F f b. 

PACS numbers: 67.85.Pq, 37.10.Jk, 11.30.Pb 



Introduction — Recently, studies in the supersymmetry 
(SUSY) for a mixture of cold Bose and Fermi atoms have 
made spectacular progress In such a cold atomic 

system, however, a Bose atom never transits to a Fermi 
atom, its superpartner or vice verse. In addition to the 
nonrelativity, this is another essential difference of this 
low-energy SUSY from the SUSY in high-energy physics. 
For the latter, such SUSY decay processes are always an- 
ticipated, e.g., a quark (lepton) may emit or absorb a 
gaugino and decays to a squark (slepton), the superpart- 
ner of the quark (lepton) [4]. 

To expose the interesting SUSY nature of the mixture, 
the effective "decay" process must be introduced. For 
a cold atomic SUSY mixture with Bose-Einstein con- 
densation, there is an effective decay of SUSY genera- 
tors since they behave as the fermion annihilation and 
creation operators [3]. Therefore, the SUSY excitations 
can be simulated by a boson-enhanced fermionic exci- 
tation. As a result, a Nambu-Goldstone-fermion-like 
(or "goldstino-like") mode in the condensation phase of 
bosons could be observed by means of the single-particle 
spectroscopy [1, |6|. 

To achieve an exact SUSY mixture, the system param- 
eters have to be fine-tuned, which requires elaborate ex- 
perimental setups and then loses the generality. In this 
article, we explore how to observe the SUSY response 
by means of a spectroscopy measurement, even if the 
mixture deviates slightly from the SUSY and the bosons 
do not condense to form a whole ordered phase. This 
can resolve the fine-tuning restraints in measuring the 
SUSY response. On the other hand, the explicit break- 
ing of the SUSY may create new excitations, the bosonic 
particle- fermionic hole individual continuous excitations, 
other than the collective goldtino-like mode. Although 
our theory is nonrelativistic, the creation of these new ex- 
citations due to SUSY explicit breaking should be quite 
general . This may be a helpful point to the study of 
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FIG. 1: (Golor online) (a) Up-panel: The optical lattice with 
cold particles. The gray (green), black (red), and white dots 
denote Fermi atoms /, F, and the molecule 6, respectively. 
Low panel: The PA processes of two atoms to one molecule 
with the binding energy Eh. (b) The Feynman diagram for the 
linear response theory: The Green function of Q is calculated 
by RPA. Up-panel: The wavy (red) lines and dotted (red) line 
denote the free Green functions of photon and F, respectively. 
Low panel: The solid (green) curves and dotted (black) curves 
denote the free Green functions of / and 6, respectively. 



SUSY in relativistic theory. 

We consider a mixture of Bose molecules h and Fermi 
atoms / with on-site interaction in a d-dimensional opti- 
cal lattice {d = 2,3) (see Fig. 1(a)). With properly tuned 
interactions and hopping amplitudes, this h-f mixture 
may become SUSY [3|. We are interested in a special 
kind of molecule 6, a bound state of / and another species 
of Fermi atom F with binding energy F5, and restrict our 
analysis to the normal phase of the h-f mixture [7]. To 
probe the SUSY behaviors, we load a free Fermi atom F 
gas, which does not interact both h and / directly. In a 
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photoassociation(PA) process [8|, the transitions between 
two atoms and one molecule, i.e., Ff ^ 6, are induced 
by two laser beams with frequencies uji and uj2- For the 
SUSY b-f mixture, this resembles a high energy physics 
process: a 'quark' or a 'lepton' (/) absorbs a 'fermionic 
gaugino' ('absorbs' a F and emits a photon) and decays 
to a 'squark' or a 'slepton' (b) or vice verse. (One can also 
consider / to be a Fermi molecule formed by the bound 
state of a Bose atom b and a Fermi atom F, i.e., processes 
Fb ^ f. We will study these processes separately.) 

For a negative detuning = CJ2 — — £^6, we show 
that the molecule dissociation process b Ff is forbid- 
den. In the formation process Ff 6, two types of new 
fermionic excitations, an individual (bosonic) particle- 
(fermionic) hole pair continuum and a collective mode, 
emerge when the SUSY in the b-f mixture is slightly bro- 
ken. For a SUSY b-f mixture, the former is completely 
suppressed while the latter in zero-monientum becomes 
an exact eigenstate, the Goldstino mode [3]. In this sense, 
we regard these excitations as the SUSY responses. The 
PA spectrum is directly related to the the molecular for- 
mation rate varying as the detuning and faithfully de- 
scribes these two types of excitations. The position of 
peak in the PA spectrum determines the frequency of 
the collective zero- momentum mode. This molecular for- 
mation rate is measured by the number variation of the 
F atoms in time. Experimentally, the number counting 
of atoms is much simpler than detecting the single atom 
spectrum. 

Model setup — The system illustrated in Fig. [TJa) is 
described by a Hamiltonian H = Hq + H^x^ where 
Ho = H})f-\-HF with H})f = H})-\-Hf-\-V . By means of the 
Feshbach resonance ^] , the scattering lengths between F 
and the b-f mixture can be adjusted to negligibly small. 
In the tight-binding approximation, one has 

i i 

where = 6^, and Fi {a = b^ f^ and F) are the 
annihilation operators of 6, / and F at site /iQ, and 
nf — are chemical potentials and the number op- 

erators at site I. The definitions of the hopping ampli- 
tudes toe and the interaction strengths Ua.^ by the Wan- 
nier function ^^^(r) can be found in the literature [iq| . 
The spatial inhomogeneity of optical lattices trapping the 
atoms and molecules has been omitted. For the subsys- 
tem b-f mixture, the Hamiltonian H^f^s SUSY invariant 
if tb = tf^ Uhh — Uhf and jih = l-if [1, 3]. In order to 
prevent the phase separation, the parameters obey, e.g., 
AntfpfUbb > U^f in 2-dimensions, where p/ is the density 
of / atoms [11]. We choose the parameters of the system 
obeying this condition. 



The PA processes are realized by simultaneously shin- 
ing two laser beams with frequencies Ui and UJ2 into the 
lattice (shown in Fig. [TJa)). The cji-beam may turn two 
free atoms / and F into a higher energy bound state |1) 
which then may transit to the molecule b by emitting 
a photon with frequency 002- Meanwhile, the molecule b 
may also be excited to |1) by the a;2-beam and then is un- 
bound with some probability by emitting a photon with 
frequency uji. For large detuning Aq, the state |1) can 
be eliminated adiabatically, so that the PA is modeled 
by the tight-binding Hamiltonian 

He. = Y^igMme'^"' + H.c), (2) 

i 

where the detuning Sq = uJc — Fb with the effective driven 
frequency c^c = ^2 — ^1, and gj = ^oexp(— iko • r^) with 
go Gc J (i^rexp(— iko • T)wl{T)wf{T)wF{'r) being the cou- 
pHng intensity independent of the site. 

In the k-space, the Hamiltonian H^x is rewritten as 

He. = 5ov^(E Qk-k„^^ke'*°* + H.c), (3) 

k 

where qI = ^p+k/p/V^ and = exp(-ik • 

rj)/\^ (where V stands for the volume and ag stand for 
^k, /k, and Fk). p = N/V is the total density of b and / 
with the particle number N = ^^{n\ + n{). 

Molecular formation rate — The formation rate of the 
molecules b can be counted by the PA variation of F- 
fermion number R = dt {i^{t) \ Np \i^{t)) for \^p{t)) being 
the time evolution from the ground state |G) = |^) |F) of 
Hq. It follows from the linear response theory that 

R = 2g^opT.^^^^{K-5o), (4) 

k 

where the retarded Green function is given by 

Dn{Kuj)= r dxA{k-ko.x) M^l'M^^) (5) 
J -00 ^ - ^k - ^ - ^0+ 

in terms of one loop calculations (Fig. [I{b)). The sin- 
gle particle dispersions are eg = —2ta^g^i cos ks ~ 
Pa where the lattice spacings are set to be the unit. 
Ufix) is the Fermi distribution at temperature T and 
the spectral function 74(k, a;) = — ImHi^(k, cj)/7r is 
defined by the retarded Green functionHi^(k, cj) = 
-i dt{g\{Q^{t),Qim\g)e'-K 

At sufficiently low temperature, the pole and branch 
cut in Eq. (|5]) are not qualitatively affected by T and nor 
is the molecular formation rate. For simplicity, we take 
a zero temperature approximation in our calculation. It 
follows from Eq. (j5j) that the rate R = Rb^Ff — Rpf^b 
contains two parts : 

Rb^Ff = ^27TglpA{k-ko,e^ - So)0{So - £^), 

k 

^Ff^b = ^27rglpA{k-ko,e^ - So)0{-e^). (6) 

k 
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which respectively are the dissociation rate for b Ff 
and the formation rate for Ff — > b. 

Collective and individual fermionic modes — In order 
to obtain R for the weak interactions, we perturbatively 
calculate Ilji{k^uj) = p^^pQ ^(k, a;) -\-U})f]~^^ which for- 
mally results from the equation of motion of Qk- It then 
follows from the random phase approximation (RPA) il- 
lustrated by the "bubble" in (Fig. mb)) that 
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where £^kp = ^k+p ~ ^p + ^Pb^U + U^fp with dU = Utb - 
Ubf and p5 = Ni)/V; ni){x) is the Bose distribution. The 
isolated pole and branch cut of ni^(k, cj) describe the 
collective and individual SUSY excitations of \g). 

Next we consider the elementary excitations in the 
two-dimensional lattice with / atoms at half filling, i.e., 
pf = 0.5. For the SUSY b-f mixture, i.e., 5U = and 
St = tb — tf = 0, the dispersion of the collective modes, 
^c(k) ^ A/i — a |k| for the small |k| [3], is read out 
from the poles of the retarded Green function IIi^(k, cj) 
(see Fig.[2fa)), where Ap = pf—pi). For large |k|, the en- 
ergy Ec(k) = £^c(|k| , ^) depends not only on |k| but also 
on the angle 6 = diTctdin{ky / kx) . The energy ^c(|k| ^6) 
of Qk 1^) decreases as |k| increases for a fixed 6. The re- 
tarded Green's function 11/^(0, cj) = {u — Ap + iO~^)~^ 
possesses a pole u = A/i, which corresponds to the 
goldstino-like excitation Qq \g). This recovers the re- 
sult in Ref. [3] . The excitation spectrum is schematically 
shown in Fig. [2fb). 

For the b-f mixture deviating slightly from SUSY, 
the retarded Green function llji{0^uj) has an isolated 
pole and a branch cut, which correspond to a collec- 
tive fermionic mode and individual (hosonic) particle- 
(fermionic) hole pair continuum modes^ respectively. 
The pole in ujq < Eq — A5t for 5t > {) (or cjq < for 
5t < 0) describes the shifted goldstino-like mode. The 
frequencies uq of collective zero-momentum mode for dif- 
ferent St are shown in Fig. Efc). For k ^ 0, the pole of 
ni^(k, cj) has the form E^^{k.) uq — |k|^ for small 
|k|. Remarkably, the branch cut Iq of 11/^(0, cj) emerges, 
which describes individual zero-momentum modes. Here, 
lo = [Eo-4.St, Eo] for St > (or [Eq, Eo-4.St] for St<0), 
and ^0 = Ap-\-2pi)SU -\-Ubfp. Notice that for the weak in- 
teractions Ubb and Ubf the SUSY breaking from SU does 
not develop a branch cut but only shifts the positions 
of the pole and branch cut. The pole and the branch 
cut can be seen in the spectral function 74(k, cj), i.e., the 
peak and the hump in Fig. Efd) for k = 0. Note that 
for the SUSY 6-/, the branch cut length /q of 11^^(0, cj) 
shrinks to zero so that the individual continuum modes 
of zero momentum are completely suppressed. On the 
other hand, as the b-f mixture deviates from the SUSY, 
the goldstino-like mode is gradually suppressed. We ex- 
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FIG. 2: (Color online) The dispersions, spectrum and spectral 
function of the excitations for Ap — 3.93 and Uhh — Uhf — 
0.1; th is taken as the unit, (a) The dispersion of collective 
mode for the SUSY mixture, (b) The spectrum for the SUSY 
system: the dashed (red) and dotted (blue) curves denote the 
dispersions of collective modes for ^ = and 7r/5, respectively. 
The solid (blue) curve denotes the dispersion of atoms F for a 
small Fermi momentum kp- (c) The frequencies of Q\ \g) for 
different St. (d) The spectral function of the zero- momentum 
for St = —0.1: The sharp peak is the shifted goldstino-like 
mode. 



amine the dependence of the spectral function on the in- 
teracting strength and find that the hump height may be 
depressed as the interaction becomes stronger, e.g., the 
height is lower than 0.5 for Uhf = Uhh = 0.5 comparing 
with - 10 in Fig. [Sfd) for Uhf = Uhh = ^-l- 

In order to study the PA spectrum of the molecular 
formation rate, we discuss the excitation spectrum shown 
in Fig. [2fb). For some momenta k, there is a collective 
mode (dashed red curve) below the individual continuum. 
For other momenta kc, the dispersion of the collective 
mode merges into the continuum. However, for small 
momentum k, there always exists a collective mode below 
the individual continuum. For convenience, we define a 
critical momentum /cq(6>), so that for a fixed when 
|k| > kqiO) the negative frequencies of the mode Qk 
emerge, i.e., 74(k, cj < 0) 7^ when |k| > kQ{0)^ and 
A(k,cj < 0) = when |k| < kgiO). 

PA spectrum — The rate R varies as detuning Sq or the 
light frequency ujc Measurement of the b boson forma- 
tion rate varying as Sq is called the PA spectrum S{So). 
For a long wave photon, the coupling gj varies slowly in 
space and the rates in Eq. (|6]) are approximately inde- 
pendent of ko. 

According to Eq. (|6]), the dissociation rate Rh^Ff does 
not vanish only if Sq > and A(k, s^ — Sq) ^ 0. Because 
the spectral function A(k, x < 0) ^ is defined by the re- 
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FIG. 3: (Color online) The PA spectra with the same pa- 
rameters as those in Fig. 2. (a) The detuning dependence 
for St = 0.05 and —0.05, the solid (red) and dashed (blue) 
lines, respectively. The unit of S is 27TgopNF and Nf/N is 
taken to be 0.01. (b) The major peak values for different St. 
This shows that the peak value is suppressed when the system 
deviates from SUSY. 



tarded Green function, it does not vanish only when the 
energies for collective modes or individual modes of Qk 
are negative for the large |k| > kg {6). We consider a di- 
lute Fermi gas F with the chemical potential /up ^ — 4^^, 
the dispersion relation turns out = ti? |k|^ — /ieff, 
where /ieff = /iF + 4ti?. In this case, the fermion F pos- 
sesses a small Fermi momentum kp = \/ j^es/tr which is 
much smaller than kQ{0) for small deviating St. There- 
fore, is always positive when |k| > kq (see Fig. [2jb)). 
For the negative detuning Sq, the condition Sq > is not 
satisfied in the regime |k| > kqiO). That is, A(k, £^ — Sq) 
and 0{So—£^) can not be non-zero simultaneously for the 
negative detuning and small Fermi momentum kp- This 
finishes our proof of R^^Ff = 0. 

The vanishing of R^^Ff for the negative Sq and small 
Fermi momentum kF can be understood in a more 
straightforward way. The transition b Ff is described 
by the Hermite conjugate term (H.c.) in the Hamilto- 
nian H^x^ which is a high-frequency oscillation term when 
^0 < 0. Hence, the Fermi golden rule results in that 
Rb^F f vanishing under the first-order perturbation (lin- 
ear response). 

For the negative Sq and small Fermi momentum kF 
(/ieff ti?), the molecule formation rate now is reduced 
to R = —RFf^b and 
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which leads to our main result: The PA spectrum 
S{So) = —RFf^b (see Fig. displays the spectral 

function of excitations qI \g). 

For the SUSY b-f mixture, the length of branch cut Iq 
tends to zero and the individual modes are suppressed. 
Meanwhile, the residue Zq = 1 and the formation rate 



RFf- 



glN/[2(tF a)] = Ro (X N at So = -A/i, 



and vanishes for the other detunings. There is a sharp 
peak at = —A/i in the PA spectrum. For a generic 
b-f mixture deviating from SUSY, the residue Zq < 1 



decreases as increases. As a result, the peak height 
is lowered while its position is shifted to Sq = —ooq. The 
ratio RFf^b{^o)/Ro is shown in Fig. [3fb) for different 
St and a small /ieff, where RFf^bi^o) is the value of 
RFf^b at So = —ujQ. Remarkably, a minor hump de- 
velops in the region Sq G Iq due to the emergence of 
individual modes (See Fig.[3fa)). As the system deviates 
further from SUSY, the individual modes are enhanced 
due to the sum rules J dujA{0^uj) = 1. The temperature 
may suppress and broaden the peak and hump. These 
characters of the PA spectrum in the b-f mixture are 
experimentally measurable SUSY responses to the light 
field. 

Conclusions — We studied how to observe the SUSY 
nature of the b-f mixture in optical lattices through PA 
spectra. For the Bose molecules formed with two species 
of Fermi atoms, we showed that the photon induced 
atom-molecule transition displays the signal of SUSY. As 
the response to the PA processes, a fermionic individual 
continuum and the goldstino-like mode were found. The 
PA spectrum can explicitly witness the molecular forma- 
tion rate of Ff b. Because the goldstino-like mode 
in zero momentum turns to be the exact eigen state for 
the SUSY mixture, the major peak in the PA spectrum 
reflects the SUSY response to the light field, even if the 
mixture is not fine-tuned to a SUSY one. 
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